      SUBROUTINE SETUP
      RETURN
      END
C   HS60
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS60'/
      DATA (X(I),I=1,3)/3*2.D0/
      DATA N/3/ , NH/1/ , NG/6/
      DATA TAU/.1D0/ , TAU0/1.D0/ , DEL0/1.D0/
      DATA (GUNIT(1,J),J=2,7)/6*1/,(GUNIT(3,I),I=2,7)/3*1,3*-1/,
     F     (GUNIT(2,I),I=2,7)/1,2,3,1,2,3/
     F     ,((GUNIT(I,J),I=1,3),J=0,1)/-1,0,0,-1,0,0/
      END
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),FX
      ICF=ICF+1
      FX=(X(1)-1.D0)**2 + (X(1)-X(2))**2 + (X(2)-X(3))**4
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),GRADF(NX)
      ICGF=ICGF+1
      GRADF(1)= 4. D0*X(1)-2.D0*X(2) - 2.D0
      GRADF(2)= 2. D0*(X(2)-X(1))+4.D0*(X(2)-X(3))**3
      GRADF(3)=-4. D0*(X(2)-X(3))**3
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
      HXI=X(1)*(1.D0+X(2)**2)+X(3)**4-4.D0-3.D0*DSQRT(2.D0)
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      GRADHI(1)=1.D0+X(2)**2
      GRADHI(2)=2.D0*X(1)*X(2)
      GRADHI(3)=4.D0*X(3)**3
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GXI
      IF ( GUNIT(1,I+NH) .EQ. -1 ) CRES(I+NH)=CRES(I+NH)+1
      IF(I.GT.3)      GXI=10.D0-X(I-3)
      IF(I.LE.3)      GXI=X(I)+1.D1
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX) ,GRADGI(NX)
       DO      50      J=1,3
       GRADGI(J)=0.D0
   50  CONTINUE
       IF(I .GT. 3)    GRADGI(I-3)=-1.D0
       IF(I .LE. 3)    GRADGI(I)=1.D0
      RETURN
      END


